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ABSTRACT 


The  characterizing  feature  of  multistatic  systems  is  the  inclusion  of 
geographically  independent  receivers  that  might  hear  echoes  from  any  source.  We 
develop  a  simple  analytic  theory  to  predict  detection  probability,  and  use  it  to  study 
cost/effectiveness  issues  and  the  advisability  of  using  colocated  source/receiver  pairs 
instead  of  independent  sources  and  receivers. 
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IV 


1. 


INTRODUCTION 


We  consider  underwater  sonar  systems  composed  of  sources  and  receivers. 
Detection  occurs  when  a  source  generates  a  brief  pulse  of  underwater  sound  that  is 
reflected  by  its  submarine  target  and  then  detected  by  one  or  more  receivers,  thereby 
revealing  the  submarine’s  location.  The  source  of  sound  might  be  a  ship,  a  helicopter 
with  a  dipping  sonar,  an  explosive  charge  dropped  by  an  aircraft,  or  an  active  sonobuoy. 
If  the  receiver  is  always  colocated  with  the  source,  the  system  is  called  monostatic,  and 
can  be  roughly  characterized  by  the  distance  at  which  detection  is  barely  possible,  the 
“detection  range”.  In  multistatic  systems,  the  sources  and  receivers  are  not  colocated. 
Multistatic  systems  are  more  difficult  to  characterize  than  monostatic  systems,  but  have 
some  arguments  in  their  favor.  Among  them  are 

•  Receivers  are  less  expensive  than  sources,  so  it  makes  sense  to  employ 
more  receivers  than  sources.  For  example,  the  U.S.  Navy’s  new  SSQ-125 
active  sonobuoy  is  expected  to  initially  cost  about  seven  times  as  much  as 
the  SSQ-1 10  passive  sonobuoy  that  will  listen  for  its  signals  (J.  Cardarelli, 
personal  communication,  February  2010). 

•  A  multistatic  system  can  employ  different  platforms  for  sources  and 
receivers.  A  ship  might  be  the  source,  while  the  receivers  are  sonobuoys. 

•  Sources  reveal  their  locations  to  submarines  when  they  transmit,  and 
submarines  can  use  that  information  to  avoid  detection.  However,  the 
independent  receivers  in  a  multistatic  system  do  not  reveal  their  positions. 

•  It  is  possible  that  the  reflected  signal  will  be  received  by  multiple 
receivers.  The  resulting  location  estimate  will  be  more  precise,  and  the 
phenomenon  can  also  help  to  eliminate  some  of  the  false  alarms  that 
monostatic  active  systems  are  normally  prey  to.  Coon  (1997)  discusses  the 
fusion  of  detections  from  multiple  receivers. 

Our  main  ambition  here  is  to  quantify  the  probability  of  detection  in  a  multistatic 
sonobuoy  system.  The  restriction  to  sonobuoys  (hereafter  simply  “buoys”)  is  because  we 
intend  to  make  no  allowance  for  the  effect  of  motion  on  the  part  of  sources,  receivers,  or 
targets — all  three  entities  are  assumed  to  remain  stationary  in  two-dimensional  space.  We 
will  include  the  possibility  that  some  of  the  sources  are  “posts”,  which  will  mean  that  the 
source  has  a  colocated  receiver.  Although  we  will  consistently  refer  to  detection  by 
sonobuoy  systems,  detection  by  multistatic  radar  systems  is  subject  to  the  same 
mathematics. 

The  effectiveness  of  any  sonobuoy  pattern  will  depend  on  the  geometric 
arrangement  of  the  buoys,  and  the  problem  of  selecting  the  arrangement  to  maximize  the 
detection  probability  therefore  arises.  We  will  not  consider  such  optimization  questions, 
assuming  instead  that  all  buoys  are  simply  located  uniformly  at  random  within  some 
region.  Our  excuse  for  this  is  that  our  goals  are  strategic  rather  than  tactical.  Specifically, 
the  goal  is  to  quantify  the  detection  probability  as  a  function  of  the  numbers  of  buoys  of 
various  types  that  are  employed,  rather  than  to  detennine  exactly  how  they  ought  to  be 
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arranged.  While  detection  probability  depends  on  the  geometric  arrangement,  it  depends 
even  more  strongly  on  resources,  and  the  assumption  of  random  deployment  will  make  it 
possible  to  expose  the  dependence  on  resources  in  an  analytically  simple  manner. 

The  U.S.  Navy  employs  ASPECT,  a  tactical  decision  aid  that  is  capable  of 
predicting  detection  probability  for  multistatic  sonobuoy  fields.  ASPECT  considers 
geometric  optimization  questions,  and  employs  a  physical  model  that  is  more  realistic 
than  anything  that  will  be  employed  here.  However,  ASPECT  is  a  menu-driven 
Monte  Carlo  simulation  that  is  not  analytically  suitable  for  considering  resource 
allocation  questions.  Bowen  and  Mitnick  (1999)  describe  the  Multistatic  Performance 
Prediction  Methodology  (MPPM).  In  MPPM,  sources  are  arranged  on  a  rectangular  grid 
with  fixed  spacing,  rather  than  randomly  distributed  as  assumed  here.  Other  multistatic 
models  include  the  Multistatic  Acoustic  Simulation  Model  (MSASM,  on  which  ASPECT 
is  based),  the  Sonar  Equation  Modeling  and  Simulation  Tool  (SE-MAST)  and  the 
Surveillance  Operational  Concepts  Model  (SOCM).  All  of  these  models  are  more 
realistic  and  less  tractable  than  the  ones  considered  here. 

Many  of  the  figures  and  tables  referred  to  below  were  generated  using  a  Microsoft 
Excel™  workbook  named  IEER.xls.  Possession  of  this  workbook  would  allow  the  reader 
to  vary  parameters,  iterate  simulations,  verify  logic,  and  possibly  discover  mistakes  made 
by  the  author,  who  would  appreciate  being  notified  if  that  is  the  case.  The  workbook  can 
be  found  among  the  downloads  at  http://faculty.nps.edu/awashburn/  . 

2.  RANDOMLY  PLACED  FIELDS  OF  INDEPENDENT  SOURCES  AND 

RECEIVERS 

2.1  Detection  Probability 

Our  model  of  “randomly  placed  field”  will  throughout  be  a  Poisson  field.  A 
Poisson  field  of  points  in  ^-dimensional  Euclidean  space  is  characterized  by  a  single 
parameter  X,  representing  the  average  number  of  points  per  unit  volume.  Poisson  fields 
have  many  desirable  analytic  properties,  one  of  which  is  that  the  number  of  points  inside 
any  region  with  volume  V  is  a  Poisson  random  variable  with  mean  XV.  In  particular,  the 
probability  of  finding  no  points  inside  such  a  region  is  exp(-A,  V),  a  fact  that  will  be 
frequently  employed  below.  We  will  also  use  two  additional  properties: 

•  The  superposition  of  a  Poisson  field  with  density  X  on  an  independent 
Poisson  field  with  density  p  is  a  Poisson  field  with  density  X  +  p. 

•  If  a  Poisson  field  with  density  X  is  “thinned”  by  removing  each  of  its 
points  independently  with  probability  p,  then  the  remaining  points 
constitute  a  Poisson  field  with  density  X(\-p). 

Consider,  then,  two  independent,  two-dimensional  Poisson  fields,  one  (sources) 
with  density  g  and  the  other  (receivers)  with  density  h.  Each  source  emits  an 
omnidirectional  sound  that  is  reflected  by  a  target,  and  the  reflected  energy  is  eventually 
received  by  each  receiver.  We  assume  that  the  transmission  loss  at  a  distance  R  is 
proportional  to  R  a  for  some  positive  value  of  a.  Spherical  spreading  corresponds  to 
a  =  2,  but  other  values  of  a  will  reasonably  model  any  direct-path  situation,  so  we  refer 
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to  this  as  a  direct-path  model.  If  Ri  and  R2  are  the  distances  of  a  source  and  a  receiver 
from  the  target,  then  the  total  transmission  loss  is  t  =  {R{a)(R{a)  =  (RlR2ya ,  so 


detection  will  depend  on  whether  the  product  RXR2  is  smaller  than  some  threshold.  The 

region  where  detection  is  possible  is  thus  the  interior  of  some  Cassinian  oval  (Cox,  1989). 
The  detection  threshold  depends  on  the  source  level,  the  target  strength,  the  sensitivity  of 
the  receiver  and  the  background  noise  level,  but  all  of  these  can  be  combined  into  a  single 
constant  p  with  dimensions  of  length  such  that  detection  happens  if  and  only  if 
RXR2  <  p1 .  Only  the  smallest  such  product  need  concern  us,  since  detection  by  any 
source-receiver  pair  is  sufficient  for  our  purposes.  We  intend  to  quantify  Q,  the 
probability  that  the  smallest  range  product  is  larger  than  p2.  Q  is  the  nondetection 
probability,  and  1  -  Q  is  the  detection  probability. 

Since  the  sources  (receivers)  fonn  a  Poisson  field  with  density  g  ( h ),  the  two  tail 
functions  of  interest  are  P(RX  >  r)  =  exp(-;r gr 2 )  and  P(R2  >  r)  =  exp(-Tz7?/-2 ) ,  for  all 

r>  0.  In  each  case,  the  event  that  the  nearest  distance  is  larger  than  r  is  the  same  as  the 
event  that  there  are  no  points  in  a  circle  with  radius  r;  that  is,  we  are  employing  the 
formula  for  the  probability  that  a  Poisson  random  variable  is  zero.  Since  R\  and  AS  are 
independent,  we  have 

oo  2  00  4 

Q  =  J P(RX  >  — )  fK  ( r)dr  =j exp(-  n^2  )2 nhr  exp {-nhr1  )dr  .  (2.1) 

o  r  o  r 


Let  s  =  7i gp2  and  t  =  nhp1  be  dimensionless  versions  of  the  source  and  receiver 
densities,  and  let  x  =  nhr1 .  Substituting  x,  5,  and  t  into  (2. 1),  we  have 

,st 
'  x 


. 

Q  =  J  exp(-(—  +  x))dx  =  yKx  (y) , 


(2.2) 


where  y  =  2\[st  and  ifi(y)  is  a  modified  Bessel  function  of  order  1  (BESSELK(y,l)  in 
Microsoft  Excel™). 

The  detection  probability  is 

P{y)  =  l-Q  =  l-yKi{yy  (2.3) 


The  detection  probability  depends  only  on  the  “effort  density”  parameter  y,  which 
incorporates  all  factors  of  tactical  relevance.  Figure  1  shows  P(y)  together  with  two 
approximations,  one  of  which  is  accurate  for  small  y  and  the  other  for  large  y.  The  two 
approximations  are 

P(y)  ~  -4rln(^)  for  small  v,  and  P(y) «  1- exP(“T)(1  +  ^-)  for  iarge  y  (2.4) 
2  2  V  2  8y 


Both  approximations  are  taken  from  Abramovitz  and  Stegun  (1964). 

Now  suppose  that  x\  sources  and  x2  receivers  are  randomly  placed  inside  a  region 
of  area  A',  with  no  buoys  outside  of  A',  and  that  A'  is  immersed  in  a  larger  region  A 
wherein  a  target  is  placed  at  random  (Figure  7).  We  wish  to  choose  A'  to  maximize  the 
probability  that  the  target  is  detected.  For  monostatic  systems,  the  best  choice  of  A'  is  to 
make  it  as  large  as  possible  (nearly  all  of  A),  since  doing  so  makes  it  unlikely  that  the 
regions  covered  by  the  individual  buoys  will  wastefully  overlap.  This  is  not  necessarily 
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true  with  multistatic  systems,  the  reason  being  that  P(y)  is  not  a  concave  function  when  y 
is  small,  as  is  evident  in  Figure  1. 

We  assume  that  the  target  will  not  be  detected  if  its  position  lies  outside  of  A',  and 
that  a  target  inside  of  A'  essentially  faces  two  infinite  Poisson  fields  with  densities 
g  =  x\!A'  and  h  =  xil  A'.  The  word  “essentially”  is  carefully  chosen,  since  the  numbers  of 
buoys  within  A'  would  be  random  in  a  Poisson  field,  rather  than  the  fixed  numbers  xi  and 
X2,  but  there  should  be  little  difference  in  effectiveness  when  those  numbers  are  large. 
The  first  assumption  is  pessimistic,  since  targets  that  are  close  to  A'  can  still  be  detected, 
and  the  second  is  optimistic,  since  targets  near  the  edge  of  A'  do  not  face  complete 
Poisson  fields.  Perhaps  the  two  assumptions  taken  together  are  neutral,  at  least  if  A'  is 
large.  This  possible  neutrality  will  be  tested  in  Section  6,  but  our  purpose  here  is  simply 
to  explore  the  consequences  of  making  both  assumptions. 
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Figure  1:  The  detection  probability  P(y),  together  with  two  approximations. 


The  probability  that  the  target  is  located  in  A'  is  A' /A,  and  within  A'  the  product  of 


which  determines  the  conditional  detection  probability.  The  searcher  can  have  a  large  y' 

in  a  small  area  or  a  small  y'  in  a  large  area,  and  must  decide  which  is  best.  The 
unconditional  detection  probability  is 


(2.5) 
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where  the  last  equality  defines  v.  The  searcher’s  problem  is  to  maximize  the  ratio  of  P(y' ) 
to  y'  by  adjusting  A'.  Figure  2  shows  a  graph  of  the  ratio.  The  maximizing  y'  is  about  1 . 1 
and  the  maximized  ratio  is  almost  exactly  0.4.  However,  y'  cannot  exceed  y,  since  A' 
cannot  exceed  A.  Accounting  for  this,  the  maximized  PD  can  be  obtained  from  (2.5). 
Letting  the  optimized  detection  probability  be  PD  *  (y) ,  we  have 


PD*{y) 


(  0.4  v;  y  <  1.1 
[P(y);y>l.r 


(2.6) 


2 

Example:  Suppose  that  p  =  1km,  and  that  the  area  to  be  searched  is  A  =  200  km  . 
Assume  that  x\  =  20  sources  and  xi  =  40  receivers  are  available,  so  that  y  =  0.89.  Then  the 
detection  probability,  according  to  (2.6),  is  PD*{ 0.89)  =  0.36,  which  is  achieved  by 
placing  all  buoys  in  a  fractional  (0.89/1.1)  part  of  A. 

If  the  number  of  receivers  is  increased  to  160,  theny  doubles  to  1.78,  so  the  buoys 
should  now  be  spread  over  all  of  A.  The  associated  detection  probability  is 
PD*(  1.78)  =  0.67. 

If  the  number  of  receivers  is  instead  decreased  to  10,  theny  shrinks  to  0.445,  and 
the  detection  probability  becomes  A7Z7>*(0.445)  =  0.18  (exactly  half  of  what  it  is  when 
there  are  40  receivers).  If  the  buoys  were  mistakenly  distributed  evenly  over  all  of  A, 
instead  of  over  the  optimal  fractional  area,  the  associated  detection  probability  would  be 
only  P(0.445)  =  0.15. 


Figure  2:  A  graph  of  P(y')/y',  the  detection  probability  per  unit  of  effort. 


2.2  Equivalent  Covered  Area 

A  pattern  of  buoys  establishes  a  probability  of  detection  for  a  target  at  point  (x,y), 
whether  that  point  is  inside  or  outside  of  the  area  within  which  the  buoys  are  deployed. 
When  integrated  over  the  whole  plane,  that  detection  function  becomes  what  might  be 
called  an  “equivalent  covered  area”  (C),  with  equivalency  in  the  sense  that,  if  a  target  is 
uniformly  distributed  over  a  region  with  area  A  that  is  much  larger  than  C,  then  the 
detection  probability  is  C/A.  The  question  of  maximizing  C  arises.  For  monostatic  buoys, 
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the  question  is  not  interesting  because  the  best  sonobuoy  field  (ignoring  practical 
questions  such  as  how  to  construct  or  monitor  widespread  fields)  would  spread  out  over 
all  of  A.  With  multistatic  buoys,  however,  the  question  is  legitimate.  The  same  argument 
as  above  leads  to  the  conclusion  that  the  effort  density  in  the  covered  part  of  an  optimal 
field  will  be  about  1.1,  and  that  the  resulting  equivalent  area  covered  will  be 

C  =  (0A)2np2y[x&.  (2.7) 

We  do  not  propose  to  do  so  here,  but  it  would  be  a  reasonable  project  to  design 
multistatic  fields  that  are  optimal  in  the  sense  of  covering  the  maximum  equivalent  area, 
but  without  our  assumption  that  the  buoys  are  scattered  randomly.  The  equivalent  area 
covered  would  of  course  be  larger  than  (2.7). 

3.  SONOBUOY  FIELDS  THAT  INCLUDE  POSTS 

There  are  some  good  arguments,  tactical  convenience  being  one  of  them,  for 
deploying  buoys  in  “posts”  that  consist  of  colocated  sources  and  receivers.  Cox  (1989,  p. 
23),  in  response  to  a  question  about  his  seminal  multistatic  work,  stated  that  “It  makes 
sense  to  me  to  use  multiple  receivers  with  a  single  high  powered  source.  One  of  these 
should  be  monostatic  so  that  TLi  =  TL2.  The  others  are  bistatic  and  provide  increased 
coverage  and  countenneasure  resistance.” 

In  principle,  one  could  have  three  types  of  buoys:  posts,  receivers,  and  sources. 
Since  we  expect  receivers  to  be  cheaper  than  sources,  however,  we  begin  by  considering 
fields  that  consist  of  posts  and  receivers,  but  no  (independent)  sources,  just  as 
Cox  envisioned. 

3.1  Sonobuoy  Fields  That  Include  Posts  and  Receivers,  but  No 
Independent  Sources 

3.1.1  Detection  Probability 

Sonobuoy  fields  of  the  type  considered  in  this  section  will  be  called  PR 
(Post/Receiver)  fields.  There  is  an  identical  theory  for  fields  that  include  independent 
sources,  but  no  independent  receivers,  since  the  roles  of  receivers  and  sources  can  always 
be  interchanged.  One  of  our  objects  is  to  compare  the  effectiveness  of  PR  fields  with  that 
of  the  SR  (Source/Receiver)  fields  that  were  considered  above  in  Section  2.  PR  fields 
have  no  independent  sources,  whereas  SR  fields  have  no  posts. 

If  any  post  in  a  PR  field  is  closer  to  the  target  than  p,  then  the  miss  probability  is 
zero.  Otherwise,  the  same  logic  that  underlies  (2.2)  applies,  so  the  formula  for  miss 
probability  in  a  PR  field  is  similar  to  (2.2).  If  s  is  the  dimensionless  density  of  posts  and  t 
is  the  dimensionless  density  of  independent  receivers,  then  the  miss  probability  is 
Q{s,st),  where 

00 

Q{s,v)  =  Jexp(-(—  +  x))dx .  (3.1) 

S  ^ 

The  analytic  effect  of  changing  sources  into  posts  is  to  change  the  lower  limit  of 
the  integral  from  0  in  (2.2)  to  s  in  (3.1),  thus  decreasing  the  miss  probability.  If  there  are 
no  independent  receivers  (t  =  0),  the  miss  probability  according  to  (3.1)  is  now  exp(-s) 
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instead  of  1.  There  seems  to  be  no  way  to  express  (3.1)  in  terms  of  commonly  available 
functions.  Even  so,  it  is  not  difficult  to  approximate  Q(s,v )  accurately  using 
numerical  integration. 

Formula  (3.1)  assumes  that  the  buoys  are  spread  out  in  a  constant  density  over  the 
whole  plane.  If  x\  posts  and  x2  independent  receivers  are  spread  out  randomly  over  a  large 
area  A,  then  we  can  estimate  the  detection  probability  within  A  by  substituting 
5  =  np2x j  /  A  and  t  =  np1  x2 1  A,  as  in  Section  2.  Also,  as  in  Section  2,  we  might  consider 

the  possibility  that  the  buoys  should  be  concentrated  in  a  smaller  area  A',  or  even  in  two 
overlapping  areas,  one  for  posts  and  one  for  independent  receivers.  However,  the 
motivation  for  such  concentration  tactics  is  less  strong  here  because  the  posts,  considered 
by  themselves,  should  theoretically  be  spread  out  over  the  entirety  of  A.  Consequently, 
we  will  not  consider  the  potential  benefits  of  buoy  concentration  in  PR  fields. 

3.1.2  Cost/Effectiveness 


We  next  consider  the  question  of  how  a  given  budget  should  be  divided  between 
sources  and  receivers,  as  well  as  the  question  of  whether  PR  fields  are  superior  to 
SR  fields. 

First,  consider  the  SR  case.  If  sources  and  receivers  cost  c\  and  c2  each,  the  cost  of 
a  buoy  field  consisting  of  x\  sources  and  x2  receivers  is  c]x]  +  c2x2 .  Since  the  detection 

probability  is  determined  by  the  product  x\X2,  it  is  a  calculus  exercise  to  conclude  that 
sources  and  receivers  should  each  consume  half  of  the  available  budget.  If  the  budget  is 


B,  the  resulting  product  is 


B 


xxx2  = 


4cjC2 


Substituting  this  into  (2.7),  we  find  that 


C  =  (0.4)— p=iL.  If  A  is  the  area  to  be  covered,  we  can  also  compute  y  =  nPj^_  ;  from 
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Ajc\c 


which  the  detection  probability  can  be  computed  via  (2.3),  or  possibly  via  (2.6)  if  buoy 
concentration  is  considered. 

In  the  special  case  where  c\  =  c2,  the  number  of  sources  in  an  SR  field  should 
equal  the  number  of  receivers.  If,  instead,  we  were  to  deploy  a  PR  field  consisting  purely 
of  posts  that  cost  2c  \  each,  then  each  post  would  cover  an  area  of  up",  and  the  total  area 

covered  would  be - .  This  exceeds  the  SR  coverage  by  20%,  so  we  have  at  least  one 

2cj 

situation  (small  budget,  equally  expensive  buoys)  where  a  PR  field  is  preferred  to  an  SR 
field.  A  PR  field  might  be  even  more  effective  if  it  included  some  independent  receivers, 
a  possibility  that  we  consider  next. 

For  subsequent  investigations,  we  impose  the  cost  constraint  s  +  ct<b.  There  is 
no  loss  of  generality  in  restricting  the  coefficient  of  s  to  be  1 .  Suppose  instead  that  each 
source  costs  c\  and  each  receiver  costs  c2,  as  before,  with  x\  and  x2  being  the  numbers  of 
sources  and  receivers  used  to  cover  an  area  of  size  A.  If  B  is  the  total  budget  for  the 
sonobuoy  field,  the  cost  constraint  would  be  cxx 1  +  c2x2  <  B ,  which  is  equivalent  to 
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xxnp 

A 


+  c. 


x2np~  Bnp 


A 


A 


If  we  now  divide  through  by  c\  and  substitute  the 


dimensionless  densities  s  and  t  for  the  ratios  involving  xi  and  X2,  we  have  a  constraint  of 

c  B  ^ 

the  proposed  form  with  c  =  —  and  b  =  — — .  Using  this  condensed  notation,  we  can 

c,  c,A 


summarize  the  optimal  balance  between  sources  and  receivers  in  an  SR  field  by 

b  b  b2  I —  b  ^ 

5  =  —,  t  =  — ,  st  =  — ,  and y  =  I'J st  = —j=  (3.2) 

2  2  c  4c  Vc 

In  a  PR  field,  if  the  cost  of  each  post  is  the  sum  of  the  costs  of  one  source  and  one 
receiver,  then  the  proper  cost  constraint  with  c  and  b  as  defined  above  is  (1  -\-c)s+ct  <  b  . 
Let  Qmin(b,c )  be  the  result  of  minimizing  Q(s,st)  subject  to  that  constraint.  For 
sufficiently  small  budgets,  the  result  of  this  minimization  will  be  that  t=  0;  that  is,  no 
independent  receivers  should  be  utilized.  To  prove  this,  let  Q\{s,t)  and  Qi{s,t)  be  the 

derivatives  of  Q(s,st )  with  respect  to  5  and  t,  respectively,  and  define 

the  function 


Wo/  dx 

~{x +  -))  —  . 

X  X 


(3.3) 


Then  Q2(s,t)  =  -H(s,t)  and  Qx(s,t)  =  -—H(s,t)-exp(-(s  +  t)) ,  as  can  be 

s 

verified  by  differentiating  Q(s,t)  with  respect  to  its  two  arguments.  For  (,v,0)  to  be 
minimizing,  it  is  necessary  that  there  exist  a  Lagrange  multiplier  X  such  that  (5,0) 
minimizes  the  expression  Q(s,st)  + l(s(\  + c)  +  tc) ,  while  simultaneously  s  =  b  l(\  +  c) . 
The  first-order  minimization  conditions  are  that  Qi(s,0)  +  /.( 1  +c)  =  0  and 
Qi{s,0)  +  Xc  >  0.  Since  Q^sA))  =  -cxp(-.s) ,  the  first  equation  requires 
X  =  exp(-s)/(l  +  c) .  Substituting  this  into  the  inequality,  we  must 
have  exp(s)H(s,  0)  <c!{\  +  c) .  Let  //(.s  )  be  the  left-hand  side  of  this  inequality: 


H(s)  =  exp(5)J  exp(-x)  —  =  exp(5)[J 


l-exp(-w) 


du  -\n(u)-y\ . 


(3.4) 


where  y  is  Euler’s  constant  (0.577  . . .).  Figure  3  shows  this  function. 

If  the  budget  is  b,  and  if  //(/)/( 1  +c))  is  smaller  than  the  cost  ratio  c/(l+c),  then  the 
entire  budget  should  be  spent  on  posts.  It  is  still  possible  that  an  SR  field  might  be  better 
than  even  the  best  PR  field,  so  in  general  we  have 

Qmvn(b,c)  =  min{exp( — -),PD*(-^)},  for  small  b  .  (3.5) 

1  +  c  dc 


Example:  Consider  the  case  ( b,c )  =  (1,1).  Since  H(0.5)  =  0.461,  which  is  smaller 
than  0.5,  the  best  PR  field  will  have  no  independent  receivers.  However,  it  is  the  second 
tenn  in  (3.5)  that  is  minimizing,  so  the  best  field  is  actually  of  type  SR.  If  b  is  reduced 
while  c  remains  1 ,  the  best  PR  field  will  still  contain  only  posts,  but  will  also  be  superior 
to  any  SR  field,  as  in  the  earlier  example  where  covered  areas  were  compared. 
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H(s) 


Figure  3:  The  function  H(s )  is  continuous,  concave,  and  increasing,  with  /T(0)=0. 

Except  for  small  budgets,  we  have  no  analytic  expression  for  Qmin(b,c),  and  our 
only  additional  observation  is  that  the  optimal  s  in  a  PR  field  will  always  be  at  least  as 

large  as  — — — .  This  is  because  that  value  maximizes  the  st  product,  and  Q(s,v)  is  a 
2(1  +  c) 

decreasing  function  of  its  first  argument. 

Example:  Figure  4  shows  a  typical  plot  of  Q(s,st )  versus  5  when  (h,c)=(  1,  0.4). 
The  graph  demonstrates  that  Qmin(l,  0.4)  =  0.43 .  The  optimal  value  of  s  is  about  0.47, 
which  exceeds  the  value  (0.36)  that  maximizes  the  st  product.  H(l/1.4)  is  about  0.5, 
which  exceeds  c/(l+c),  so  it  should  not  be  surprising  that  independent  receivers  are 
utilized  in  this  example,  as  well  as  posts. 


Figure  4:  PR  miss  probability  as  a  function  of  s  when  b= 1  and  c=0.4. 
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Figure  5  shows  the  miss  probability  for  the  best  SR  and  PR  fields  when  all  buoys 
are  equally  expensive,  with  PR  winning  for  small  budgets.  The  SR  curve  is  based  on  (2.6) 
with  y  =  b  (y  and  b  are  equal  when  c=  1),  while  the  PR  curve  is  based  on  (3.1).  The 
initial  slope  of  the  SR  curve  is  -0.4.  The  initial  part  of  the  PR  curve  employs  no 
receivers,  and  is  consequently  equal  to  exp(-b/(l+c)),  the  initial  slope  of  which  is  -0.5. 
This  comparison  of  initial  slopes  explains  why  PR  is  superior  when  the  budget  is  small 
and  all  buoys  are  equally  expensive.  The  initial  slope  of  the  SR  curve  is,  in  general, 
-0.4/ -Jc ,  while  the  initial  slope  of  the  PR  curve  is  -l/(l+c).  By  comparing  the  two 
expressions,  it  is  not  difficult  to  show  that  the  PR  field  will  be  superior  for  small  budgets 
as  long  as  0.25  <c<  4.  For  more  extreme  values  of  c,  the  SR  field  is  superior  for  all 
budgets.  The  SR  field  is  superior  for  large  budgets  even  when  c  =  1,  as  Figure  5  makes 
clear.  Depending  on  circumstances,  either  type  of  field  can  be  superior.  If  the  SSQ-125 
active  buoy  turns  out  to  be  as  expensive  as  estimated  in  the  introduction  (c  =  1/7),  then, 
regardless  of  budget,  an  SR  field  will  be  more  cost/effective  than  a  PR  field. 

Lack  of  universal  dominance  of  either  the  PR  or  the  SR  type  suggests  that  the 
optimal  sonobuoy  field  might  actually  include  all  three  of  the  fundamental  buoy  types: 
posts,  free  receivers,  and  free  sources.  This  possibility  is  not  covered  by  any  of  the  above 
developments,  so  we  examine  it  in  the  next  section. 


Figure  5:  Miss  probability  for  SR  and  PR  fields  when  sources  and  receivers  have  the 

same  cost  (c=l). 

3.2  Sonobuoy  Fields  That  Include  Posts,  Independent  Receivers,  and 
Independent  Sources 

3.2.1  Detection  Probability 

In  the  general  case,  the  density  of  posts  is/  the  density  of  free  sources  is  g,  and 
the  density  of  free  receivers  is  h.  Let  R/,  Rg,  and  R/,  be  the  distances  from  the  target  to  the 
nearest  post,  free  source,  and  free  receiver,  respectively.  The  nondetection  event  can  be 
decomposed  into  two  mutually  exclusive  parts  E\  and  E2,  where 
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Ex  =  (Rg  <  p), (Rf  >  p), (RfRg  >  p2),{RgRh  >  P2URfRh  >  P 2) 
e2=(R>  p),(Rf>  p),(RfRe>  P2),(ReRh>  P2URfRh>  P2) 


(3.6) 


The  commas  separating  events  in  (3.6)  should  be  read  “and”.  The  target  will  not  be 
detected  if  and  only  if  there  is  no  post  within  p  of  the  target,  and  if  every  product  of  two 
ranges  exceeds  p".  These  are  the  last  four  events  in  both  E\  and  E2,  and  the  first  event 
simply  partitions  on  whether  Rg  is  small  (E\)  or  large  (£3)  (nondetection  is  also  possible  if 
Rg  =  p,  but  that  possibility  is  safely  ignored  in  (3.6)  because  it  has  probability  zero). 

In  the  expression  for  E\  in  (3.6),  the  second  event  can  be  omitted  because  it  is 
implied  by  the  first  and  third,  and  the  fifth  event  can  be  omitted  because  it  is  implied  by 
the  first,  third  and  fourth.  In  E2,  the  third  event  can  be  omitted  because  it  is  implied  by  the 
first  and  second.  Now  let  Rs  and  R,  be  the  distances  from  the  target  to  the  nearest  source 
(free  or  not)  and  receiver  (free  or  not).  Thus,  Rs  =  min(£; ,  /fjand  Rr  =min (Rf,Rh). 

Random  variables  Rs  and  R,  are  easily  described  mathematically,  since  the  posts,  together 
with  either  the  independent  sources  or  independent  receivers,  constitute  Poisson  fields 
with  densities  f+g  (sources)  and  f  +  h  (receivers).  Furthermore,  given  the  above 
observations,  E\  and  E2  can  be  expressed  as 


Ex  =  (Rg  <  p),(RrRg  >  P2)  and  E2  =  (Rs  >  p),(RsRh  >p2). 
After  conditioning  on  Rg  in  E\  and  Rs  in  E2,  we  have 

P  4 

P(Ef  =  J 2^xgexp(-^gx2)exp(-^-(/  +  h)^-f)dx,  and 

A  ^ 


P(E2)  =  J 2 7rx(f  +  g ) exp {-n(J  +  g)x2) exp(-7rh ^r)dx  . 


(3.7) 


(3.8) 


(3.9) 


2  2  2 

Now,  let  5  =  rfp  ,  t  =  7Tgp~,  and  z  =  7thp~  be  the  dimensionless  densities  of  posts, 

a 

sources,  and  receivers,  respectively,  and  substitute  u  =  7igx  into  (3.8).  We 
then  have 


P(EX)  =  J  exp {-{u  +  S  +  Z))du  . 


Similarly,  substituting  u  =  rtf  +  g)x2  into  (3.9),  we  have 
r  z(s  +  t ) 

P(E2 )  =  [  exp(-(w  +  — - ))du  =  Q(s  +  t,z(s  +  t)) , 

J  11 


(3.10) 


(3.11) 


where  Q()  is  the  function  defined  in  (3.1).  Define  the  sum  of  (3.10)  and  (3.11)  to  be 
Q(s,  t,  z): 

S+Z)¥„+fexp(-(„  +  z(s  +  ') 
u 


s  +  z 

Q(s,  t,  z)  =  J  exp(-(i/  +  - ))du  +  J  exp(-(i/  + : 


-))du . 


(3.12) 


0  "  s+t 

This  function  represents  the  miss  probability  when  a  target  is  subjected  to 
multistatic  detections  by  posts,  sources,  and  receivers  in  the  dimensionless  densities  5,  t, 
and  z,  respectively.  The  arguments  t  and  z  can  be  interchanged,  although  the  form  of 
(3.12)  does  not  make  that  apparent. 
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3.2.2  Optimization 


While  we  have  no  proof  of  the  fact,  it  appears  from  extensive  experimentation 
that  the  best  sonobuoy  field  never  includes  all  three  types,  whatever  the  budget  or  the 
costs  of  sources  and  receivers.  Given  x\  sources  and  xi  receivers,  together  with  the  ability 
to  deploy  any  pair  as  a  post,  the  best  sonobuoy  field  will  always  have  as  many  posts  as 
possible  (a  PR  field)  or  none  at  all.  Except  in  establishing  this  conclusion,  formula  (3.12) 
appears  to  have  little  use. 

4.  RELIABILITY  AND  DIRECT  BLAST 

The  idea  that  detection  is  merely  a  matter  of  being  close  enough  to  the  target 
underlies  all  of  the  above  calculations,  but  is  only  approximately  true  in  the  real  world. 
One  reason  for  this  is  that  some  buoys  may  not  function  correctly,  in  which  case  their 
distances  from  the  target  are  immaterial.  This  possibility  is  easily  handled  theoretically  in 
SR  fields.  If  rs  and  rr  are  the  reliabilities  of  sources  and  receivers,  it  is  simply  a  matter  of 
replacing  the  dimensionless  densities  5  and  t  by  srs  and  tty,  respectively.  This  is  because 
Poisson  fields  “thinned”  in  this  manner  remains  Poisson  fields.  A  similar  argument  will 
not  work  for  PR  fields  because  a  post  might  get  effectively  turned  into  an  independent 
source  if  its  receiver  should  not  function. 

Proximity  might  not  be  sufficient  even  when  all  buoys  function  as  intended.  In 
multistatic  sonar  systems,  receivers  hear  the  transmitted  signal  directly,  in  addition  to  the 
signal  reflected  from  the  target.  This  “direct  blast”  is  actually  necessary  for  locating  the 
target  because  the  difference  in  time  between  the  direct  and  reflected  signals  establishes 
an  ellipse  upon  which  the  target  must  he.  However,  the  direct  blast  is  much  stronger  than 
the  reflected  signal,  and  may  completely  obscure  it  if  the  time  difference  in  arrivals  is 
small  enough.  The  principal  situation  where  this  happens  is  when  the  target  is  more  or 
less  directly  between  the  source  and  the  receiver.  A  simple  approximation  of  the  effect  is 
the  “dead  zone”  shown  in  Figure  6,  where  receivers  will  not  detect  the  target  because  of 
the  direct  blast  arriving  nearly  simultaneously.  The  angle  9  measures  the  extent  of  the 
pie-shaped  dead  zone,  with  9  =  n  corresponding  to  no  dead  zone  at  all. 

It  is  not  difficult  to  include  the  effect  of  dead  zones  in  a  Monte  Carlo  simulation 
that  begins  by  simulating  the  locations  of  the  target,  the  sources  and  the  receivers.  For 
each  source/receiver  pair  that  passes  the  proximity  test,  one  simply  checks  whether  the 
receiver  is  in  the  dead  zone  of  the  source,  declaring  the  detection  attempt  to  be  a  failure  if 
so.  The  resulting  simulation  is  only  slightly  more  complicated  than  one  without  dead 
zones.  Unfortunately,  however,  this  simplicity  does  not  extend  to  analytic  attempts  to 
find  a  generalization  of  (2.3).  We  are  not  aware  of  any  exact  formulas  for  detection 
probability,  but  the  following  theorems  at  least  offer  lower  bounds. 

Theorem  1:  In  an  SR  field,  let  s  and  t  be  the  dimensionless  densities  of  sources 
and  receivers,  respectively,  as  in  Section  2,  but  add  the  requirement  that  the  receiver  buoy 
in  a  successful  pair  must  not  be  in  the  dead  zone  of  its  source.  Fet  p  =  6  /  n ,  and 

let  v  =  2yjpst .  If  v  is  substituted  into  (2.3)  one  has  a  lower  bound  on  the  detection 
probability. 
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Proof:  Let  E  be  the  event  that  the  target  is  detected  by  the  source  that  is  closest  to 
it.  E  will  happen  if  and  only  if  there  is  at  least  one  eligible  receiver  that  is  sufficiently 
close  to  the  target;  a  receiver  being  “eligible”  if  and  only  if  it  is  not  in  the  dead  zone  of 
the  closest  source.  The  probability  that  any  receiver  is  eligible  is  p  for  each  independent 
receiver,  so  eligible  receivers  are  a  Poisson  field  with  dimensionless  density  pt.  To 
calculate  P{E ),  we  can  now  proceed  as  in  the  derivation  of  (2.3),  except  that  pt  needs  to 
substituted  for  t.  The  effect  of  this  is  that  y  is  modified  as  in  the  statement  of  the  theorem, 
and  P(E)  =  1  -  yK^  (y) .  This  is  a  lower  bound  on  the  detection  probability  because  it  is 
possible  for  the  target  to  be  detected  even  when  the  closest  source  fails.  QED 


Figure  6:  Illustrating  a  dead  zone  where  receivers  are  useless  on  account  of 

direct  blast. 

Theorem  2:  In  a  PR  field,  let  s  and  t  be  the  dimensionless  densities  of  posts  and 
receivers,  respectively,  as  in  Section  3,  but  add  the  requirement  that  the  receiver  buoy  in  a 
successful  pair  must  not  be  in  the  dead  zone  of  its  source.  Let  p  =  0  /  n .  Then 
1  -  Q(s,  pst )  is  a  lower  bound  on  detection  probability,  where  Q()  is  the  function  defined 
in  (3.1). 

Proof:  If  the  post  that  is  nearest  to  the  target  is  closer  than  p,  then  detection  is 
certain  because  a  post’s  own  receiver  cannot  be  in  a  dead  zone.  Let  E  be  the  event  that 
the  nearest  post  is  farther  away  than  p,  but  nonetheless  detects  the  target  through  some 
other  eligible  independent  receiver.  E  will  happen  if  and  only  if  there  is  at  least  one  such 
receiver  that  is  sufficiently  close  to  the  target.  The  probability  that  any  independent 
receiver  is  eligible  is  p,  so  eligible  receivers  are  a  Poisson  field  with  dimensionless 
density  pt.  To  calculate  P(E),  we  can  now  proceed  as  in  the  derivation  of  (3.1),  except 
that  pt  needs  to  substituted  for  t.  The  result  of  this  is  as  stated  in  the  theorem.  This  is  a 
lower  bound  on  the  detection  probability  because  it  is  possible  for  the  target  to  be 
detected  even  when  the  closest  post  fails.  QED 

In  both  theorems,  dead  zones  are  handled  as  if  the  independent  receivers  had  a 
reliability  of  p.  Receivers  that  are  part  of  a  post  are  in  effect  assumed  to  be  perfectly 
reliable  in  Theorem  2. 

The  question  arises  as  to  whether  the  upper  bounds  defined  in  Theorems  1  and  2 
are  sharp.  We  note  that  the  bound  in  both  theorems  is  exact  when  9  =  n  or  (trivially) 
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when  0  =  0.  For  intermediate  values,  the  sharpness  question  will  be  tested  by  simulation 
in  Section  6.2. 

5.  CAVEATS 

We  emphasize  that  all  of  the  above  analysis  applies  only  to  the  direct-path 
transmission  mode.  The  primary  sonar  modes  where  such  a  transmission  loss  model  does 
not  apply  are  the  bottom-bounce  and  convergent-zone  (CZ)  modes.  Consider  the  latter 
case,  and  assume  that  there  is  only  one  CZ  at  a  distance  R*  (typically  about  50  km)  with 
width  A  (typically  about  5  km).  There  will  be  a  CZ  detection  if  and  only  if  the  circular 
annulus  about  the  target  at  distance  R*  and  width  A  contains  at  least  one  buoy  of  each 
type.  Since  the  area  of  this  annulus  is  (approximately)  2nAR*,  the  probability  of  this  is 

PD  =  (1  -  exp(-2;zAR  *  g))(l  -  exp(-2;zAR  *  h))  ~  (2;zAR*)2  gh  .  (5.1) 

The  approximation  in  (5.1)  depends  only  on  the  product  gh,  so  seemingly  there  is 
an  analogy  with  the  direct-path  case.  However,  the  approximation  is  actually  an  upper 
bound  that  is  accurate  only  when  PD  is  small.  For  sonobuoy  fields  that  are  dense  enough 
to  assure  a  large  detection  probability,  the  CZ  case  differs  significantly  from  the  direct 
path  case. 

Example:  Suppose  that  R*= 50  km,  A=5  km,  g=  0.001/km2,  and  h  =  0.002/km2. 
Then,  using  (5.1),  PD  =  0.76.  Suppose  there  is  also  a  direct  path  mode  with  p=  10  km. 
The  associated  y  is  about  0.89,  and  the  associated  direct  path  detection  probability  from 
(2.3)  is  about  0.36,  significantly  smaller  than  the  CZ  probability.  When  a  CZ  exists,  it 
will  generally  provide  most  of  the  detection  probability,  as  it  does  here.  However,  CZs 
can  be  ephemeral  enough  to  make  tracking  difficult  even  in  monostatic  systems,  and 
dependence  on  multiple  CZs  in  the  multistatic  case  is  unlikely  to  improve  the  situation. 

There  are  other  assumptions  made  above  that  can  differ  significantly  from  the  real 
world.  Direct  blast  regions  only  approximately  resemble  the  conical  region  of  Figure  6. 
The  pervasive  assumptions  of  radial  symmetry  made  above  can  be  significantly  false  in 
the  ocean,  where  transmission  loss  often  depends  on  direction,  as  well  as  range.  Our 
analysis  makes  no  allowance  for  target  motion,  or  for  any  phenomenon  that  might  make 
p  a  function  of  time.  Better  forecasts  of  detection  probability  in  a  particular  instance  will 
always  come  from  simulations  that  account  for  some  of  these  phenomena. 

6.  CALIBRATION  FOR  SR  FIELDS 

6.1  The  Neutrality  Assumption 

In  this  section,  we  test  the  neutrality  assumption  lying  behind  (2.5)  and  (2.7). 
Recall  that  A'  is  the  area  within  which  the  buoys  are  placed,  while  A  is  the  potentially 
larger  area  within  which  targets  are  to  be  found  (Figure  7).  The  neutrality  assumption  is 
that  targets  detected  outside  of  A'  will  compensate  for  targets  not  detected  inside  of  A'. 
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The  assumption  is  clearly  optimistic  if  A  and  A'  are  equal,  since  in  that  case  there  are  no 
targets  outside  of  A'.  It  is  just  as  clearly  pessimistic  if  buoys  are  very  dense,  since  in  that 
case  all  targets  outside  of  A'  that  are  within  p  of  A'  will  be  detected,  in  addition  to  all 
targets  inside  of  A'.  Much  depends  on  the  density  of  buoys  and  the  relationship  between 
A  and  A'.  We  begin  by  assuming  that  A  is  much  larger  than  A',  and  that  A'  has  been 
adjusted  to  maximize  the  detection  probability;  that  is,  we  wish  to  compare  the  theoretical 
equivalent  covered  area  C  from  (2.7)  with  an  experimental  version  Cexp. 

To  estimate  Cexp,  without  loss  of  generality  we  first  take  p  =  I .  We  also  take 
A'  =  2k^xxx1  / 1 . 1 ,  since  this  is  the  area  that  theoretically  maximizes  detection 

probability.  All  buoys  are  placed  within  a  square  of  side  L  with  that  area,  so  L~=A' .  Since 
detection  is  impossible  for  targets  more  than  1  unit  away  from  this  square,  we  generate  30 
targets  at  random  within  a  square  of  side  L+2  (so  A  =  (L+ 2)2),  and  count  the  number  X 
that  are  detected.  A  target  is  detected  if  the  product  of  its  distances  from  the  closest 
receiver  and  the  closest  source  is  smaller  than  1.  This  experiment  is  repeated  3,000  times; 
that  is,  all  the  buoys  are  randomly  relocated  in  A'  and  each  such  pattern  is  tested  against 
30  targets  randomly  located  in  A.  Let  X  be  the  total  number  of  targets  detected  in  all 
3,000  replications.  Then  the  equivalent  covered  area  is  Cexp  =  AX  / 90,000 . 


: .vi ; sources ; and; as  receivers 


are  placed  at  random  in  the: 
inner  square >1' ; : ; : ; : ; : ; : ; : ; : ;  ; : ; : ; 


the  target  is  placed  at  random  in  the  outer  square  A 

Figure  7:  A  coordinate  framework  for  simulation  experiments. 


Table  1  shows  the  ratio  Cexp/C  for  various  values  of  X|  and  X2.  It  should  be  evident 
that  the  neutrality  assumption  is  optimistic,  but  only  by  a  small  amount.  The  average  of 
all  100  ratios  is  96.6%,  so  the  theoretical  area  coverage  is  about  3%  high.  Some  idea  of 
the  variability  of  results  can  be  obtained  by  comparing  symmetric  entries  (xi^)  and 
(x2,xi),  since  both  have  the  same  theoretical  mean. 


15 


5  10  15  20  25  30  35  40  45  50  AVG 

0.959 
0.963 
0.962 
0.965 
0.969 
0.968 
0.967 
0.969 
0.972 
0.971 
0.966 

Table  1:  Rows  and  columns  are  labeled  by  the  number  of  sources  and  receivers, 
respectively.  The  table  entries  are  the  ratio  Cexp/C,  the  ratio  of  experimental  to 

theoretical  equivalent  covered  area. 

We  next  consider  the  case  where  A  and  A'  are  both  the  same  square,  in  which  case 
we  expect  the  theory  to  be  optimistic.  Without  loss  of  generality,  we  take  A'  to  be  a  unit 
square,  and  place  all  buoys  and  the  target  within  it.  The  fundamental  experiment  is  now 
to  measure  Z,  the  product  of  the  distances  from  the  target  to  the  nearest  source  and  the 
nearest  receiver,  recording  a  detection  if  Z  <  p1 .  Alternatively,  let  Y  =  2;r-s/x1x2  Z  and 

y  =  2n^xpc2  p 1 ,  so  that  detection  is  equivalent  to  the  event  Y  <  y.  If  we  now  measure  the 

cumulative  distribution  function  of  7,  we  can  deal  with  all  values  of  p  simultaneously. 
The  fundamental  experiment  is  repeated  with  30  independent  targets  for  each  buoy 
configuration,  and  then  for  3,000  random  buoy  configurations.  The  results  are  shown  in 
Figure  8  for  20  sources  and  40  receivers,  where  “data”  (the  sample  cumulative 
distribution  function  or  sample  CDF)  is  compared  with  “theory”  (formula  (2.3)  for  P(y)). 
It  can  be  seen  that  the  theory  is  indeed  optimistic  in  this  case.  The  contrast  would  be  even 
larger  if  the  sample  CDF  were  compared  with  (2.6)  instead  of  (2.3),  since  theory  would 
advise  concentrating  the  buoys  in  a  smaller  area  when  y  is  small.  Results  for  other 
numbers  of  sources  and  receivers  are  visually  identical,  so  are  not  shown. 

Figure  9  shows  a  similar  comparison  for  a  case  where  the  targets  are  located  in  an 
areaZ  that  is  larger  thanZ'.  Specifically,  A  has  side  of  1.2  whileZ'  still  has  a  side  of  1,  so 
A' I  A  is  1/(1. 2)  ,  or  0.694.  This  ratio  is  an  upper  bound  on  the  theoretical  detection 
probability  because  the  theoretical  assumption  is  that  targets  outside  of  A'  will  not  be 
detected.  If  p  >  0. 1 ,  however,  every  point  in  A  (except  in  the  corners)  is  within  p  of  some 
point  in  A',  so  a  sufficiently  large  density  of  buoys  within  A'  will  detect  nearly  every 
target.  Figure  9  shows  that  the  sample  CDF  agrees  closely  with  the  theoretical  formula 
(2.3)  for  small  values  of  the  effort  density  within  A',  but  that  (2.3)  predicts  too  small  a 
detection  probability  in  dense  fields  (large  values  ofy).  With  20  sources,  40  receivers, 
and  p  ~  0.1,  y'  is  1.78 — approximately  the  value  at  which  the  two  curves  begin 
to  diverge. 


5 

0.974 

0.957 

0.966 

0.960 

0.952 

0.963 

0.955 

0.955 

0.957 

0.951 

10 

0.946 

0.976 

0.959 

0.972 

0.969 

0.961 

0.955 

0.965 

0.958 

0.965 

15 

0.957 

0.957 

0.958 

0.972 

0.956 

0.968 

0.962 

0.967 

0.955 

0.966 

20 

0.963 

0.953 

0.968 

0.965 

0.972 

0.963 

0.972 

0.963 

0.965 

0.962 

25 

0.956 

0.963 

0.965 

0.968 

0.976 

0.975 

0.968 

0.967 

0.968 

0.983 

30 

0.958 

0.964 

0.963 

0.971 

0.967 

0.963 

0.976 

0.964 

0.977 

0.971 

35 

0.947 

0.966 

0.960 

0.967 

0.966 

0.971 

0.973 

0.969 

0.978 

0.977 

40 

0.960 

0.969 

0.971 

0.958 

0.972 

0.974 

0.968 

0.972 

0.976 

0.969 
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0.962 

0.965 

0.967 

0.971 

0.970 

0.975 

0.985 

0.970 

0.974 

0.976 
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0.952 

0.972 

0.967 

0.978 

0.974 

0.975 

0.969 

0.971 

0.973 

0.978 
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0.964 

0.965 

0.968 

0.967 

0.969 

0.968 

0.966 

0.968 

0.970 
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Figure  8:  Comparing  theoretical  and  experimental  detection  probabilities  in  a  case 
where  20  sources  and  40  receivers  are  randomly  distributed  over  the  same  area  that 

contains  the  target. 


Figure  9:  Comparing  theoretical  and  experimental  detection  probabilities  in 
a  case  where  20  sources  and  40  receivers  are  randomly  distributed  in  an  area  ,4'  that 
is  smaller  than  the  area  that  contains  the  target.  With  a  high  effort  density  in  A', 
some  targets  outside  of  A'  are  detected  experimentally,  but  not  in  theory. 
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6.2  A  Direct  Blast  Experiment 

Each  number  in  Table  2  is  the  ratio  of  an  experimentally  detennined  effective 
area  of  coverage  to  the  theoretical  number  determined  from  (2.7),  as  in  Table  1.  The 
experiment  is  as  in  Section  6.1,  except  that 

•  pxi  is  substituted  for  x2  in  calculating  the  effective  area  covered  from 
(2.7),  as  well  as  in  calculating  A',  where p  =  d/n; 

•  in  each  replication,  the  sonobuoy  field  is  generated  independently  within 
A'  for  each  target;  and 

•  an  angular  test  is  added  to  the  proximity  test  before  a  detection  is 
recorded,  as  pictured  in  Figure  6. 

Ten  thousand  replications  are  carried  out  in  a  Monte  Carlo  simulation  where  each 
(source,  receiver)  pair  of  buoys  is  tested  until  some  pair  finally  detects  the  target,  in 
which  case  the  detection  counter  is  incremented  and  a  new  replication  begun,  or  no 
detection  is  recorded  if  no  pair  detects  the  target. 

The  case  9  =  n  is  covered  by  Table  1,  which  shows  that  the  theoretical  formula 
(2.7)  is  about  3%  optimistic.  When  0  =  3ti/4,  Table  2  shows  that  the  formula  is  about  5% 
pessimistic.  When  9  =  n/2,  Table  3  shows  that  the  formula  is  about  13%  pessimistic,  and 
when  9  =  rc/4,  Table  4  shows  that  the  fonnula  is  about  27%  pessimistic.  Presumably  the 
increasing  pessimism  is  because  sources  other  than  the  closest  source  become 
increasingly  important  as  the  direct  blast  effect  becomes  stronger. 

AVG 
1.04 
1.05 
1.04 
1.04 
1.04 
1.05 
1.04 
1.04 
1.04 
1.04 
1.04 

Table  2:  Rows  and  columns  are  labeled  by  the  number  of  sources  and  receivers, 
respectively.  The  table  entries  are  the  ratio  Cexp/C  based  on  10,000  repetitions  of  the 
experiment  of  locating  a  target  in  an  SR  field  of  20  sources  and  40  receivers, 

with  0  =  3n/4. 


5.00  10.00  15.00  20.00  25.00  30.00  35.00  40.00  45.00  50.00 


5.00 

1.06 

1.05 

1.03 

1.04 

1.01 

1.03 

1.04 

1.02 

1.03 

1.04 

10.00 

1.04 

1.05 

1.05 

1.03 

1.04 

1.07 

1.07 

1.04 

1.05 

1.04 

15.00 

1.01 

1.07 

1.02 

1.03 

1.04 

1.04 

1.06 

1.05 

1.03 

1.04 

20.00 

1.04 

1.04 

1.02 

1.05 

1.06 

1.02 

1.04 

1.05 

1.07 

1.05 

25.00 

1.05 

1.03 

1.05 

1.03 

1.06 

1.05 

1.05 

1.02 

1.02 

1.05 

30.00 

1.03 

1.06 

1.05 

1.04 

1.03 

1.05 

1.06 

1.07 

1.06 

1.04 

35.00 

1.02 

1.06 

1.04 

1.07 

1.04 

1.05 

1.05 

1.02 

1.04 

1.05 

40.00 

1.01 

1.05 

1.03 

1.03 

1.02 

1.01 

1.04 

1.05 

1.07 

1.08 

45.00 

1.00 

1.05 

1.03 

1.04 

1.03 

1.06 

1.06 

1.05 

1.03 

1.04 

50.00 

1.02 

1.01 

1.05 

1.04 

1.05 

1.06 

1.02 

1.05 

1.06 

1.03 

AVG  1.03  1.05  1.04  1.04  1.04  1.04  1.05  1.04  1.05  1.05 
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5.00  10.00  15.00  20.00  25.00  30.00  35.00  40.00  45.00  50.00  AVG 

1.12 

1.13 

1.13 

1.13 

1.13 

1.14 

1.13 

1.13 

1.13 

1.13 

1.13 


5.00  10.00  15.00  20.00  25.00  30.00  35.00  40.00  45.00  50.00  AVG 

1.23 

1.25 

1.27 

1.27 

1.27 

1.27 

1.27 

1.28 

1.26 

1.27 

AVG  1.24  1.26  1.26  1.27  1.27  1.27  1.27  1.27  1.27  1.26  1.26 


5.00 

1.29 

1.23 

1.25 

1.22 

1.21 

1.24 

1.21 

1.21 

1.18 

1.20 

10.00 

1.26 

1.23 

1.25 

1.27 

1.29 

1.26 

1.25 

1.26 

1.24 

1.22 

15.00 

1.26 

1.25 

1.28 

1.28 

1.27 

1.29 

1.28 

1.28 

1.25 

1.28 

20.00 

1.21 

1.28 

1.28 

1.29 

1.29 

1.27 

1.27 

1.29 

1.29 

1.25 

25.00 

1.24 

1.29 

1.26 

1.26 

1.28 

1.28 

1.28 

1.25 

1.27 

1.27 

30.00 

1.22 

1.26 

1.25 

1.27 

1.29 

1.29 

1.29 

1.30 

1.30 

1.26 

35.00 

1.20 

1.26 

1.29 

1.28 

1.24 

1.29 

1.29 

1.28 

1.29 

1.26 

40.00 

1.26 

1.29 

1.24 

1.30 

1.28 

1.31 

1.28 

1.25 

1.30 

1.28 

45.00 

1.24 

1.28 

1.25 

1.26 

1.25 

1.26 

1.25 

1.28 

1.28 

1.27 

50.00 

1.25 

1.25 

1.28 

1.27 

1.29 

1.27 

1.28 

1.30 

1.29 

1.28 

5.00 

1.13 

1.14 

1.10 

1.15 

1.10 

1.12 

1.14 

1.10 

1.14 

1.12 

10.00 

1.14 

1.14 

1.18 

1.12 

1.11 

1.12 

1.12 

1.12 

1.12 

1.12 

15.00 

1.12 

1.09 

1.15 

1.12 

1.16 

1.12 

1.15 

1.13 

1.14 

1.15 

20.00 

1.15 

1.13 

1.13 

1.12 

1.11 

1.13 

1.13 

1.13 

1.13 

1.09 

25.00 

1.10 

1.12 

1.11 

1.14 

1.12 

1.13 

1.16 

1.14 

1.13 

1.15 

30.00 

1.09 

1.13 

1.14 

1.13 

1.18 

1.14 

1.15 

1.14 

1.15 

1.14 

35.00 

1.09 

1.11 

1.18 

1.16 

1.12 

1.12 

1.13 

1.11 

1.15 

1.16 

40.00 

1.10 

1.11 

1.13 

1.12 

1.14 

1.14 

1.16 

1.15 

1.15 

1.14 

45.00 

1.11 

1.09 

1.14 

1.14 

1.13 

1.14 

1.15 

1.13 

1.16 

1.13 

50.00 

1.08 

1.10 

1.11 

1.13 

1.14 

1.17 

1.13 

1.14 

1.16 

1.12 

AVG 

1.11 

1.12 

1.14 

1.13 

1.13 

1.13 

1.14 

1.13 

1.14 

1.13 

Table  3:  As  in  Table  2,  but  with  0  =  7t/2. 


Table  4:  As  in  Table  2,  but  with  0  =  n/4. 


7.  SUMMARY 

We  have  developed  several  simple  formulas  for  approximating  the  detection 
probability  of  a  multistatic  sonobuoy  field,  all  of  which  depend  on  the  monostatic 
detection  distance  p.  Either  an  SR  field  (no  posts)  or  a  PR  field  (all  sources  are  posts) 
may  be  the  most  cost/effective,  depending  on  the  budget  level  and  the  relative  costs  of  the 
two  buoy  types.  If  x\  sources  and  xi  receivers  are  used  to  detect  a  target  within  an  area  A, 
then  the  dimensionless  coverage  ratio, 

2np2  Jx.Xi 
A 

is  fundamental.  If  y  is  small  compared  to  1,  then  detection  within  A  is  unlikely,  or  the 
opposite  is  true  if  the  ratio  is  large. 
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